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Abstract 


Two  problems  related  to  the  Royal  Thai  Navy's  recruitment  effort  are 
modeled  and  solved  by  mathematical  programming.  The  first  problem  is  the 
optimal  assignment  of  draftees  to  branch  naval  bases,  which  is  formulated 
as  a  transportation  problem.  The  second  problem  is  the  optimal  routing 
of  ships  for  transporting  draftees  from  the  branch  bases  to  the  main 
base,  which  is  solved  with  an  exact  integer  programming  formulation. 


Introduction 


There  is  compulsory  military  service  for  men  in  Thailand.  Inductions 
take  place  four  times  a  year,  with  the  men  who  live  in  the  coastal  provinces 
assigned  to  the  Navy. 

When  a  man  is  called  into  the  Thai  Navy,  he  first  reports  to  a  drafting 
center  in  his  home  locality.  He  is  then  transported  by  land  vehicle  to  a 
naval  base.  If  he  is  from  a  northern  province,  this  base  is  the  main  naval 
base  in  Sataheep,  near  Bangkok.  If  he  is  from  a  southern  province,  the  draftee 
is  first  brought  overland  to  a  branch  naval  base  and  then  transported  by  ship 
to  the  main  base.  (See  the  map  of  Figure  1.)  In  the  southern  provinces  there 
are  36  drafting  centers  and  4  branch  bases.  This  paper  addresses  two  problems 
relating  to  the  assignment  and  transportation  of  draftees  from  the  southern 
provinces: 

Problem  1:   How  many  men  from  each  drafting  center  should  be  assigned  and 

transported  to  each  branch  base? 

Problem  2 :  Given  an  available  fleet  of  ships  stationed  at  the  main  base, 

which  of  those  ships  should  be  used  and  how  should  they  be  routed  so  as 

to  pick  up  all  the  men  at  the  branch  bases  and  transport  them  to  the  main 

base? 

The  second  problem  is  the  more  difficult. 

We  shall  present  a  linear  programming  model  which  is  used  for  solving  the 
first  problem  and  an  integer  programming  model  which  is  used  for  solving  the 
second.  The  two  problems  are  not  independent  because  the  solution  to  Problem 
1  determines  the  transport  requirements  for  Problem  2. 
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\  would  be  possible  to  build  a  unified  model  to  solve  both  problems  simultane- 
ously, this  would  be  inappropriate  under  the  Thai  Navy's  current  organiza- 
tional structure.  The  personnel  decisions  of  Problem  1  and  the  ship  routing 
decisions  of  Problem  2  are  made  by  different  people  in  different  places  at 
different  times. 

Problem  1:   Assigning  Draftees  to  Branch  Bases 

Problem  1  is  modeled  as  a  standard  transportation  problem. 
We  are  given: 
bj  =  the  number  of  men  to  be  transported  from  drafting  center  j, 
a±   =  the  capacity  (number  of  men)  of  branch  base  i, 
cij  =  cos*  per  roan  for  transport  to  base  i  from  center  j, 
and  we  must  determine 

Xij  =  the  number  of  men  to  be  transported  to  base  i  from  center  j. 
The  model  is 


min  S  cijxij 


1J 

s.  t. 

J 

ai, 

all 

i 

IXij     = 

i 

bJ> 

all 

j 

with  all  x-n  >  o.   The  vehicles  used  for  transporting  the  men  are  sufficiently 
small  so  that  it  is  not  necessary  to  use  a  more  elaborate  model  (e.g,  a  fixed- 
charge  transportation  model)  for  Problem  1.   The  results  of  Problem  1  which 
are  inputs  to  Problem  2  are 

di  =  £x:Lj 
i 

0 

which  is  the  number  of  men  assigned  to  branch  base  i.   The  data  and  numerical 

results  for  a  specific  instance  of  Problem  1  (Summer,  1983)  are  given  in 

Tables  1  and  2. 
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TABLE  1.   DATA  FOR  TRANSPORTATION  PROBLEM 


Drafting 

Number  of  Men 

Center 

Distances 

to  Branch  Bases 

at  Center 

Bl 

B2 

B3 

B4 

PC 

25 

100 

170 

270 

98 

KR 

30 

85 

165 

260 

82 

PN 

10 

87 

150 

250 

123 

S\V 

23 

7^ 

140 

240 

54 

RN 

50 

75 

150 

240 

118 

CY 

77 

20 

90 

180 

98 

TK 

112 

29 

87 

175 

73 

TP 

130 

70 

120 

200 

66 

BD 

100 

0 

75 

170 

102 

SC 

125 

50 

40 

140 

57 

TS 

130 

60 

23 

125 

61 

PG 

160 

7~ 

112 

175 

104 

KL 

170 

90 

120 

178 

60 

TL 

150 

110 

120 

170 

58 

PK 

195 

120 

125 

70 

98 

KB 

175 

82 

80 

138 

107 

CW 

150 

50 

28 

123 

92 

TG 

170 

"5 

25 

100 

107 

RB 

175 

8C 

20 

87 

93 

KM 

180 

80 

48 

85 

62 

RT 

210 

125 

73 

5u 

69 

PL 

215 

125 

70 

50 

88 

TR 

210 

120 

^5 

77 

93 

KT 

220 

127 

82 

78 

So 

SN 

227 

135 

74 

28 

74 

HT 

267 

175 

110 

12 

103 

PR 

260 

167 

112 

70 

52 

ST 

282 

185 

132 

55 

77 

TA 

282 

208 

133 

37 

59 

SD 

280 

190 

130 

31 

68 

PT 

302 

220 

150 

60 

83 

YR 

310 

228 

162 

73 

61 

SB 

325 

235 

175 

81 

82 

NW 

348 

270 

200 

107 

96 

YL 

325 

230 

168 

70 

72 

3T 

360 

275 

203 

100 

53 

Base 

Capacities 

1000 

1100 

1300 

1500 

TABLE  2. 

OPTIMAL  SOLUTION  FOR 

TRANSPORTATION  PROBLEM 

From 

To 

Number  of  Men 

Cost /Man 

Extension 

PC 

Bl 

98 

25 

2450 

KR 

£1 

82 

30 

2460 

PN 

Bl 

123 

10 

1230 

SW 

Bl 

54 

23 

1242 

RN 

bl 

118 

50 

5900 

CY 

B2 

98 

20 

1960 

TK 

32 

73 

29 

2117 

TP 

z2 

66 

70 

4620 

BD 

B2 

102 

0 

0 

PG 

B2 

104 

77 

8008 

KL 

3., 

60 

90 

5400 

TL 

B2 

58 

110 

6380 

PK 

B2 

98 

120 

11760 

SC 

B3 

57 

40 

2280 

TS 

b3 

61 

23 

1403 

KB 

b3 

107 

80 

8560 

OW 

33 

92 

28 

2576 

TG 

33 

107 

25 

2675 

Re 

B3 

93 

20 

1860 

KM 

B3 

62 

48 

2976 

TR 

B3 

93 

"5 

6975 

RT 

B4 

69 

50 

3450 

PL 

S~ 

88 

50 

4400 

KT 

i34 

86 

78 

6708 

SN 

5~ 

74 

28 

2072 

HY 

B4 

103 

12 

1236 

PR 

B4 

52 

70 

3640 

ST 

B4 

77 

55 

4235 

TA 

3^ 

59 

37 

2183 

SD 

B4 

68 

'- 

2108 

PT 

B4 

83 

60 

4980 

YR 

B4 

61 

n3 

4453 

SB 

B4 

82 

81 

6642 

NW 

34 

96 

107 

10272 

YL 

B4 

72 

70 

5040 

BT 

B4 

53 

100 

5300 

Total  Cost  =  149,551 


Base 


Number  of  Men  Assigned 


Bl 

475 

32 

659 

B3 

672 

B4 

1123 

Problem  2:  Routing  Ships  for  Transporting  Draftees  to  the  Main  Base 

The  problem  of  optimally  deploying  the  available  fleet  of  ships  to  transport 
draftees  to  the  main  base  is  a  type  of  vehicle  routing  problem.  There  is  a  large 
number  of  available  approaches  to  these  problems.  (See,  e.g.,  Eell  et.  al. 
[19831  and  Bodin  et.  al.  [19831).  The  model  given  below  is  related  to  the  set- 
partitioning  approach,  in  that  we  first  create  a  list  of  possible  routes  that 
ships  can  take  and  then  model  the  problem  in  terms  of  variables  representing  the 
number  of  times  each  route  is  used.  Since  ships  are  not  identical,  these  var- 
iables must  also  denote  which  type  of  ship  is  chosen. 

In  most  vehicle  routing  problems,  the  objective  function  is  the  total  dis- 
tance traveled  or  some  concomitant  cost  measure.  The  Thai  Navy's  situation  is 
different.  Their  primary  objective  is  to  minimize  the  number  of  ships  required 
for  transporting  all  the  draftees.  This  is  because  defending  the  country,  not 
transporting  draftees,  is  the  principal  mission  of  the  Navy  fleet.  As  few  ves- 
sels as  possible  should  be  assigned  to  this  duty,  so  that  the  remaining  ships  can 
be  available  for  defense.  Since  we  anticipate  that  numerous  solutions  would  al- 
ternatively optimize  this  primary  objective,  we  use  distance  minimization  as  a 
secondary,  tie-breaking  objective. 

The  fleet  available  for  draftee  transport  at  the  time  of  this  study  con- 
sisted of  9  vessels  of  3  different  classes.  This  fleet  is  described  in 
Table  3- 

TABLE  3.   FLEET  DESCRIPTION 

Class  Capacity  of  Ship  Number  of         Bases  that  Class 

of  Ship         (Number  of  Men)  Ships  in  Class        Can  Visit 

1  100  2  1 

2  200  3  1,  2,  3 

3  600  4  1,  2,  3,  4 


From  this  information  we  generated  23  possible  voyages  for  the  model  to 

consider  using.  A  voyage  is  defined  as  a  particular  class  of  ship  assigned  to 

sail  a  particular  route.   The  23  voyages  are  described  in  Table  4.   Node  0  in 
the  route  specifications  means  the  main  base. 

TABLE  4.   VOYAGE  SPECIFICATION 
Voyage       Class        Route  Distance  (km)     Capacity  (men) 


1 

1 

0-1-0 

370 

100 

2 

2 

0-1-0 

370 

200 

3 

2 

0-2-0 

460 

200 

4 

2 

0-3-0 

600 

200 

5 

2 

0-1-2-0 

515 

200 

6 

2 

0-1-3-0 

640 

200 

7 

2 

0-2-3-0 

665 

200 

8 

2 

0-1-2-3- 

■0 

720 

200 

9 

3 

0-1-0 

370 

600 

10 

3 

0-2-0 

460 

600 

11 

3 

0-3-0 

600 

600 

12 

3 

0-4-0 

750 

600 

13 

3 

0-1-2-0 

515 

600 

14 

3 

0-1-3-0 

640 

600 

15 

3 

0-1-4-0 

810 

600 

16 

3 

0-2-3-0 

665 

600 

17 

3 

0-2-4-0 

305 

600 

18 

3 

0-3-4-0 

800 

600 

19 

3 

0-1-2-3- 

•0 

720 

600 

20 

3 

0-1-2-4- 

■0 

860 

600 

21 

3 

0-1-3-4- 

>J 

340 

600 

22 

3 

0-2-3-4- 

•0 

865 

600 

23 

3 

0-1-2-3- 

■4-0 

920 

600 

The  method  used  for  generating  the  set  of  available  voyages  is   the  fol- 
lowing.  For  each  ship  class,  we  enumerate  every   subset  of  bases  that  can  be 
visited  on  a  voyage  and  then  we  compute  the  shortest  route  for  that  subset. 
All  voyages  worth  considering  are  thereby  generated.   If  the  ship  routing  prob- 
lem were  posed  by  another  fleet  whose  number  of  classes  and  number  of  bases 
were  much  greater  than  the  Thai's,   then  this  exact  approach  might  become  com- 
putationally infeasible. 


The  variables  in  our  integer  programming  model  for  Problem  2  are: 

zj  =  number  of  times  voyage  j  is  used,  and 

y-H  =  number  of  men  transported  from  base  i  to  the  main  base  via 
voyage  j. 

The  data  for  the  integer  programming  model  are: 
Vj[  =  set  of  voyages  that  visit  base  i, 
S]£  =  set  of  voyages  that  use  a  ship  of  class  k, 
Bj  =  set  of  bases  visited  on  voyage  j, 
dj_  =  number  of  men  to  be  transported  from  base  i, 
n^  =  number  of  ships  available  of  class  k, 
cj  =  capacity  of  voyage  j,  and 
rj  =  distance  traveled  by  voyage  j  in  tens  of  thousands  of  kilometers. 

The  model  is 

min  E  (1  +  r .  )  z  . 
j      J   J 
s.  t. 

.£„   y. .  =  di,   all  i 

.Ex,   y.  .  <  c  .z  . ,  all  j 

ieBj  °ij    j  j'     J 

.  E0   z .   <  nv,   all  k 
z .  >   0  and  integer 

The  objective  function,  as  noted,  primarily  serves  to  minimize  the  number 
of  voyages  required.  The  rj  have  small  values  and  serve  only  to  break  ties. 
The  first  constraint  set  insures  that  all  the  drafted  men  at  base  i  are  trans- 
ported, using  voyages  that  in  fact  pass  trough  base  i.  The  second  constraint 
set  insures  that  the  number  of  men  transported  via  voyages  of 


type  j  is  within  the  capacity   that  will  be  made  available.   The  third  con- 
straint set  insures  that  each  ship  in  the  fleet  is  used  at  most  once.   A  ship 
can  not  be  reused  because  all  the  draftees  have  to  start  their  training  at  the 
main  base  at  the  same  time. 

Software  Employed 
The  transportation  model  of  Problem  1  is  of  course  a  standard  problem  for 
which  numerous  effective  packages  exist.  In  our  case  we  used  GNET,  the  ca- 
pacitated transshipment  problem  solver  of  Bradley,  Brown  and  Graves  [1977],  We 
wrote  a  problem  generator  and  report  writer  in  FORTRAN  to  use  in  conjunction 
with  this  program,  and  we  used  the  DEC-10  computer's  TOPS-10  operating  system 
for  managing  the  necessary  file  interfaces. 

The  integer  programming  model  in  the  particular  instance  we  solved  was 
small  enough  to  handle  with  the  general  purpose  solver  UNDO  of  Schrage 
[1981].  The  model  formulation  above  for  the  Thai  fleet  of  Tables  3  and  4  has 
23  integer  variables  Zi  and  45  continuous  variables  yi j  •  The  zj  variables  can 
take  on  any  integer  value  up  to  and  including  the  number  of  ships  in  the  class 
assigned  to  voyage  j.  LINDO,  like  many  other  commercial  integer  programming 
solvers,  requires  that  all  the  integer  variables  be  binary.  Therefore,  we  had 
to  reformulate  the  model  with  the  zj  replaced  by  the  binary  expansions: 

z  .  =  z  .,  +  2z  .0     •  =  1     # 

Zj  =  Zjl  +  2zJ2  +  4zj3,   j  =  9,.. .,23 

This  resulted  in  a  model  with  61  binary  variables. 

Another  problem  generator  was  written  to  automatically  generate  the 
voyages  and  create  this  input,  so  that  other  instances  of  the  model  can  be 
treated  routinely.   The  results  for  the  model  instance  based  on  Tabls  3  and 
4  are  given  in  Tables  5  and  6. 


TABLE  5.   OPTIMAL  SOLUTION  TO  SHIP  ROUTING  MODEL 
Nonzero  Variable        Value  Interpretation 


Z(2,2) 

1 

Z(5,l) 

1 

ZdO,l) 

1 

Z(ll,l) 

1 

Z(12,l) 

I 

Z(18,l) 

1 

Y(l, 2) 

400 

Yd, 5) 

75 

Y(2,5) 

59 

Y(2,10) 

600 

Y(3,1D 

595 

Y(3,18) 

77 

Y(4,12) 

600 

Y(4,18) 

523 

Use  voyage  2  twice. 

Use  voyage  5  once. 

Use  voyage  10  once. 

Use  voyage  11  once. 

Use  voyage  12  once. 

Use  voyage  18  once. 

Transport  200  men  on  each  voyage  2. 

Transport  75  men  from  Bl  on  voyage  5. 

Transport  59  men  from  B2  on  voyage  5. 

Transport  600  men  from  B2  on  voyage  10. 

Transport  595  men  from  B3  on  voyage  11. 

Transport  77  men  from  B3  on  voyage  18. 

Transport  600  men  from  B4  on  voyage  12 . 

Transport  523  men  from  B4  on  voyage  18. 


Class 

Voyage 

of  Ship 

Route 

Men 

Carried 

2 

2 

0-1-0 

200 

2 

2 

0-1-0 

200 

5 

: 

0-1-2-0 

75+59 

10 

3 

0-2-0 

600 

11 

3 

0-3-0 

595 

12 

3 

0-4-0 

600 

18 

3 

0-3-4-0 

77+523 

TABLE  6.   OPTIMAL  SHIPPING  SCHEDULE 

Distance 

370 
370 
515 
460 
600 
750 
800 

Total  number  of  voyages  =  7 

Total  distance         =  3865  km 

Discussion  of  Results 
According  to  Thai  Naval  authorities,  the  results  obtained  by  the  math- 
ematical programming  models  are  superior  to  the  results  obtained  through  the 
manual  procedures  currently  employed.  In  the  instance  of  the  transpor-tation 
model  reported  above,  it  turned  out  that  just  sending  each  man  to  the  nearest 
base,  as  would  be  done  by   the  manual  procedure,   turned  out   to  be 
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feasible.  However,  in  other  instances,  this  procedure  has  resulted  in  vio- 
lations of  the  base  capacities,  which  necessitated  additional  transport ion  of 
men  from  oversubscribed  bases  to  undersubscribed  bases.  The  authorities  were 
pleased  to  see  that  the  cost  and  time  delays  of  these  extra  trips  could  be 
avoided. 

With  respect  to  the  ship  routing  model,  again  the  authorities  were  pleas- 
ed to  see  our  recommendation  of  an  efficient  deployment  of  their  fleet.  Data 
gathering  is  currently  in  process  to  determine  how  many  voyages  and  kilometers 
could  have  been  saved  if  the  mathematical  programming  models  had  been  avail- 
able for  the  past  few  years. 

Conclusion 

Two  mathematical  programming  models,  a  standard  transportation  model  and 
a  new  and  exact  integer  programming  formulation,  have  been  shown  to  be  ef- 
fective for  optimally  routing  ships  and  assigning  personnel  for  Thai  Naval 
recruitment.  The  models  were  solved  by  general  purpose  software.  Our  ap- 
proach is  potentially  applicable  to  other  situations,  however  it  is  likely 
that  special  purpose  software  would  be  needed  for  the  integer  programming  pro- 
blem if  larger  fleets  are  considered. 
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